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on  the  propagation  of  an  ultrasonic  wave:  it  decreases  the  stiffness  and  increases  the 
attenuation.  The  central  objective  of  this  work  has  been  to  correlate  damage  states  with 
changes  in  the  two  ultrasonic  parameters  (wavespoed  and  attenuation). 

We  have  developed  a  new  technique  for  measuring  the  wavespeed  and  attenuation  in  the 
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the  excitation  and  detection  of  lamo  waves  in  the  lengthwise  direction.  Thus,  both  the 
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the  stiffness  was  found  to  be  insensitive  to  transverse  cracking.  The  attenuation, 
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K- 

It  Is  well-known  that  composite  materials  develop  a  complex  damage  state 
when  they  are  subjected  to  monotonic  or  fatigue  loading.  The  damage  has.  In 
general,  two  effects  on  the  propagation  of  an  ultrasonic  wave:  it  decreases 
the  stiffness  and  increases  the  attenuation.  The  central  objective  of  this 
work  has  been  to  correlate  damage  states  with  changes  in  the  two  ultrasonic 
parameters  (wavespeed  and  attenuation). 

We  have  developed  a  new  technique  f o**  measuring  the  wavespeed  end 
attenuation  In  the  thickness  direction.  In  extremely  thin  laminates.  We  have 
also  developed  a  technique  for  the  excitation  and  detection  of  Lamb  waves  in 
the  lengthwise  direction.  Thus  both  the  in-plane  and  out-of-plane 
measurements  can  be  made. 

Damage  in  the  form  of  transverse  cracking  In  cross-ply,  laminates  has 
been  studied  by  the  use  of  these  two  techniques.  For  through-the-thlckness 
measurements  the  stiffness  was  found  to  be  insensitive  to  transverse 
cracking.  The  attenuation,  however,  was  found  to  be  quite  sensitive  and, 
therefore,  has  been  shown  to  be  a  reliable  damage  metric.  For  the 
complementary  case  of  Lamb  wave  propagation  in  the  lengthwise  direction,  both 
the  stiffness  and  the  attenuation  were  observed  to  be  sensitive  damage 
parameters. 


1.  SUMMARY  OF  RESEARCH  ACCOMPLISHED 

The  following  Is  a  brief  suninary  of  the  research  performed  during  the 
grant  period: 

1.1  A  theoretical  Investigation  of  plane  wave  propagation  In  a  general 
unbounded  anisotropic  solid. 

1.2  A  combined  theoret icdi /experimental  study  of  plane  wave  propagation.  In 
the  thickness  direction,  in  a  composite  plate  immersed  In  water. 

1.3  A  combined  theoretical/experimental  study  of  Lamb  wave  propagation,  in 
the  lengthwise  direction.  In  a  composite  plate  immersed  In  water. 

1.4  A  study  of  damage  in  the  form  of  transverse  cracking  by  the  use  of  these 
two  techniques 
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2.  DETAILS  OF  RESEARCH  ACCOMPLISHED 
2.1  Wave  Propagation  In  an  Anisotropic  Solid 

As  a  prelude  to  the  design  of  experimental  techniques  and  the 
interpretation  of  the  gathered  data  It  was  considered  desirable  to  study,  in 
fairly  general  terms,  the  nature  of  wave  propagation  in  anisotropic  media.  In 
particular,  several  wave  propagation  concepts  which  are  unique  to  anisotropic 
media  were  studied  in  detail:  velocity  surface,  wave  surface,  slowness 
surface,  group  velocity,  and  deviations  of  the  energy  flux  vector  and  particle 
displacement  vector  from  the  wave  vector.  These  results  were  suimiarized  in  an 
interim  progress  report  to  the  AFOSR  [1].  In  the  following  we  reproduce  a 
typical  result  of  [11.  A  generic  unidirectional  f iber-reinforced 
graphite/epoxy  composite  is  considered;  the  fibers  are  along  the  x^  axis 
(transverse  isotropy  with  x^  as  the  axis  of  Isotropy).  Fig  1  shows  the 
velocity  surface  i.e.  the  distance  between  the  origin  and  any  point  on  the 
curve  Is  the  velocity  of  a  plane  wave  In  that  particular  direction.  The 
velocity  is  an  eigenvalue  of  the  well-known  Christoffel  equations 

^Ci jkinjni  ‘  °v^6 1k^  pk  *  0  W 

where  Is  the  stiffness  matrix,  o  is  the  density,  pk  is  the  unit  particle 
displacement  vector  (eigenvector)  and  v  Is  the  phase  velocity  of  time-harmonic 
waves  travelling  In  the  direction  n^.  Eq  (1)  has  three  eigenvalues:  one 
quasi-longitudinal  (L)  and  two  quasi-shear  (SH  and  5V)  .  Note  that,  as 
expected,  the  highest  wavespeed  is  that  of  the  longitudinal  wave  in  the  fiber 
direction.  The  unit  vector  p  is  thc  particle  displacement  vector.  This 
brings  out  clearly  a  unique  feature  of  wave  propagation  in  anisotropic 
materials:  particle  displacement  in  general  is  not  parallel  (perpendicular)  to 
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the  wavevector  for  longitudinal  (transverse)  waves.  In  turn,  this  gives  rise 
to  another  Interesting  phenomenon,  namely,  the  direction  in  which  the  energy 
flows  Is  not  in  the  same  as  the  direction  of  wave  propagation.  This  is  shown 
In  Fig.  2.  Let  slowness  s  *  1/v.  The  slowness  surfaces  for  the  longitudinal 
wave  and  two  shear  waves  are  shown  in  Fig.  2.  It  is  well-known  that  the 
energy  flux  vector  (F)  Is  perpendicular  to  the  slowness  surface;  we  have 
plotted  F  as  a  bi-directional  vector  every  ten  degrees  on  the  slowness 
surfaces.  With  reference  to  the  L-surface  we  note  that  when  the  wave  is 
propagating  anywhere  In  the  range  0°  -  80°  from  the  fiber  direction  the  energy 
flows  (practically)  in  the  fiber  direction.  This  was  perhaps  the  most 
significant  observation  of  our  theoretical  investigation. 

We  have  developed  a  computer  code  which  can  compute  various  three- 
dimensional  surfaces  for  the  most  general  (21  elastic  constants)  composite 
material. 

The  concept  of  stress-wave-factor  originated  by  Alex  Vary  [21  has 
recently  been  applied  to  composite  materials  [31 .  To  the  best  of  our 
knowledge  there  has  not  appeared  a  satisfactory  explanation  as  to  why  the  SWF 
technique  works.  Our  work  has  offered  some  (albeit  far  from  complete) 
explanation  of  this  phenomenon  [1|. 

2.2  Throuqh-The-Thickness  Measurements.  Theory 

Our  objective  here  Is  to  measure  the  wavespeed  c  and  the  attenuation  k£ 
of  a  laminate  in  its  thickness  direction.  A  conventional  method  is  the 
toneburst  method  [4].  However  the  toneburst  method  is  suitable  only  for 
relatively  thick  specimens  (about  five-wavelength  thick).  The  central 
objective  of  the  development  of  composite  materials  is  to  fabricate  light¬ 
weight  aerospace  structures  i.e.  to  use  thin  laminates  (of  the  order  of  a  few 
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millimeters);  here  the  toneburst  method  completely  breaks  down.  This 
observation  motivated  us  to  define  the  first  goal  of  the  present  research:  to 
develop  a  new  experimental  technique  suitable  for  thin  plates  (sub-millimeter 
in  thickness). 

We  have  used  the  following  fact  from  the  theory  of  Fourier  Transforms:  If 
two  events  are  very  close  In  the  time  domain,  then  a  pair  of  corresponding 
events  are  proportionately  far  apart  in  the  frequency  domain. 

In  order  to  make  our  technique  useful  in  other  areas  of  mechanics  and 
physics  of  solids,  throughout  this  investigation  we  have  adapted  a  "black-box" 
approach  l.e.  one  can  measure  the  phase  velocity,  the  group  velocity  and  the 
attenuation  of  any  linear  viscoelastic  solid.  The  technique  works  equally 
well  for  thick  or  thin  plates,  for  dispersive  or  non-disperslve  materials. 
The  details  of  this  technique  have  been  given  In  15).  Consider  a  linear 
viscoelastic  plate  (e.g.  a  composite  laminate)  Immersed  in  an  elastic  fluid 
(e.g.  water).  A  plane-fronted  longitudinal  wave  begins  to  travel  at  time  t=0 
In  the  positive  x  direction;  see  Fig.  3  for  a  Lagranglan  diagram.  Let  the 
displacement  in  the  incident  field,  Ray  1,  be  given  by 

uinc  »  f0(wt-k0x) 

where  w  is  the  circular  frequency  and  kQ  Is  the  wavenumber  in  water;  the  speed 
of  sound  in  water  cQ  =  w/kQ.  The  Interaction  of  this  pulse  with  the  plate 
results  In  a  series  of  reflected  pulses,  Rays  2,  6,  10,  ...®,  and  a  series  of 
transmitted  pulses,  Rays  4,  8,  12 . .  see  Ref  !5)  for  details. 

We  recall  the  Fourier  transfers  of  a  function  f(t)  as 

F*(u»)  =  (2*)'14  J  f(t)  e"1u,tdt,  — <«<- 
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Fig.  3  lagrangian  diagram  for  reflected  and  transmitted  field 
from  a  plate  immersed  in  water. 


with 


— co<£  <a> 


f(t)  =  (2t)_l5  jV(«)  etiwtd«. 


2.2.1  Thick  Specimens 

The  use  of  the  words  “thick"  and  "thin"  is,  of  course,  quite  arbitrary. 

For  the  purpose  of  this  discussion  we  define  a  specimen  to  be  thick  if,  in 

Fig.  3,  we  can  clearly  separate  various  pulses  from  each  other.  The  Incident 

wave  Is  the  signal  produced  by  a  typical  ultrasonic  transducer  when  energized 

by  a  short-duration  electrical  spike. 

Let  us  first  consider  the  reflected  field.  Let  f(t)  be  the  signal 

corresponding  to  Ray  2  and  g(t)  be  the  signal  corresponding  to  Rays  2  and  6 
★  ★ 

combined;  let  F  (<d),  G  (w)  be  their  Fourier  transforms,  then  (5) 


G  (uj)  _  ,  T  t  _ — 1 2k  h 
777 ' 1  •  Ti2T2ie 

F  (id) 


where  T^2  Is  the  transmission  coefficient  from  water  Into  specimen,  T2^  Is  the 
transmission  coefficient  from  specimen  into  water,  h  is  the  specimen  thickness 
and  k  Is  the  complex-valued  wavenumber  in  the  specimen,  k  *  kj  i-  1  k2,  where 
the  phase  velocity  c  *  w/kj  and  k2  Is  the  attenuation  (nepers/mm)  l.e.  a  wave 

_  U  v 

decays  exponentially  as  e  2  .  The  stiffness,  In  turn.  Is  given  by  an 

2 

expression  of  the  type  E  =  pc  .  Thus,  the  measurement  of  k  Is  the  central 

*  it 

objective  of  this  work.  If  one  writes  -(G  /F  1n  1ts  P°1ar  form  as 

Me1*,  then  It  can  be  readily  shown  that 


k1  «  -*/2h 
k2  =  (In  i)/2h 
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Furthermore,  a  plot  of  jG*(w)|  vs  u  is  characterized  by  a  series  of  resonance 
peaks  whose  spacing  Is  given  by 

Aw  =  (Cn)/h 


As  was  mentioned  earlier  this  is  a  fundamental  observation  of  this 
work.  As  the  plate  becomes  very  thin  (h*0)  the  front  and  back  surface 
reflections  become  very  close  in  the  time  domain;  however.  Aw*®  l.e.  any  two 
neighboring  resonance  peaks  in  the  frequen'.y  domain  (the  corresponding  pair  of 
events)  become  very  far  apart. 

Let  us  now  consider  the  transmitted  field  for  a  thick  specimen.  Let  f(t) 
be  the  Incident  field  with  the  specimen  removed.  Let  g(t)  be  the  signal  due 
to  Ray  4  alone.  Then 


* 

—  =  T 
F 


T  e"i(kh  +  koh) 
12  21e  0 


2.2.2  Thin  Specimens 

Consider  a  case  when  h  Is  so  small  that  Rays  2,  6,  10  ...  cannot  be 
separated  in  the  time  domain.  As  previously  mentioned,  none  of  the  existing 
techniques  can  be  used  to  measure  the  phase  velocity  or  attenuation  for  this 
case. 

Let  g(t)  be  the  total  reflected  field,  Rays  2+6+10+  ...®  .  Let  f(t)  be 
Ray  2  obtained  separately  (see  I  SI  for  details).  Then  the  complex  wavenumber 
k  can  be  obtained  from  the  set  of  equations 


-  12kh 
e 


2. 
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Z  =  P  (1+S). 
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ft, 2^21 

1 12 *21 


1) 


Wa  now  consider  the  transmitted  field.  Let  f(t)  be  the  incident  field 
through  water  only  i.e.  with  the  specimen  removed.  Let  g(t)  be  the  total 
transmitted  field,  Rays  4+8+12+..-  '.  Unlike  the  cases  considered  so  far,  here 
we  get  a  quadratic  in  Z  =  exp(-ikh) 


ZZ  +  ZY  -  Dq  =  0,  where 


Y 


1 12T2 1  F 
R21Zo  G 


★ 

★ 


ZQ  *  exp(-1hko) 


The  fact  that  one  gets  two  roots  for  the  wavenumber  is  not  at  all 
disturbing.  For  the  correct  root  the  phase  of  z  decreases  with  frequency  (for 
the  Incorrect  root  the  converse  is  true) 

This  sums  up  the  development  of  the  theory  in  support  of  the  experimental 
work  which  Is  described  next. 
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2.3  Through-the-Thlckness  Measurement.  Experimental 

A  schematic  of  the  apparatus  Is  shown  In  Fig.  4.  The  heart  of  the  system 
is  a  pair  of  accurately  matched,  broad-band,  water-immersion  piezoelectric 
transducers.  The  frequency  was  varied  in  the  range  0.5  -  1C  MHz;  the 
oscilloscope  digitizes  at  the  rate  of  100  x  10^  points  per  second.  We 
estimated  the  precision  of  our  experiments  as  follows:  For  homogeneous 
materials  (e.g.  aluminum,  epoxy,  steel,  plexiglas,  etc.)  the  precision  In 
velocity  Is  about  0.1*.  For  graphite/epoxy  composite  the  precision  in 
velocity  is  0.2*  and  In  attenuation  it  is  about  5*. 

Calibration  Procedures  Whenever  a  new  technique  for  measuring  wavespeed  (or 
attenuation)  Is  developed,  properly  it  should  be  used  to  measure  wavespeed  of 
a  material  whose  wavespeed  Is  known  to  an  accuracy  about  ten  times  better  than 
the  (claimed)  accuracy  of  the  new  experimental  technique.  Unfortunately, 
however,  we  found  that  the  National  Bureau  of  Standards  has  not  yet  developed 
a  reference  standard  for  wavespeed  (or  stiffness).  This  section  Is  devoted  to 
our  efforts  In  lieu  of  a  standard  calibration. 

ASTM  Round  Robin  Test 

Our  laboratory  participated  in  a  six- laboratory  round-robin  program 
conducted  by  the  American  Society  for  Testing  and  Materials.  Two  nickel-based 
alloys  were  tested.  We  measured  both  the  longitudinal  as  well  as  the  shear 
wavespeeds.  We  make  two  observations. 

1.  The  precision  in  our  measurements  was  ±0.15* 

2.  Measurements  by  all  six  laboratories  were  within  ±2.5*  of  each  other 

The  results  of  this  Investigation  have  been  accepted  for  publication  by  the 
Journal  of  Testing  and  Evaluation  ( 6 1 . 
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Fig.  4  Block  diagram  cf  the  apparatus. 
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Since  the  goal  of  our  effort  Is  to  develop  a  technique  for  testing  “thin" 
plates,  we  subjected  our  procedures  to  the  following  critical  test.  We 
started  with  a  "thick"  aluminum  plate  (2.807  mm)  and  gradually  machined  It 
down  to  an  extremely  "thin"  plate  of  0.258  mm  thickness  (10x10“^  Inches  i.e.  a 
very  thin  foil).  At  each  stage  we  measured  the  wavespeed.  The  results  are 
presented  in  Table  1.  (Here,  o  Is  the  standard  deviation  of  ten  measurements 
far  each  thickness).  We  conclude  that  as  we  transition  from  "thick"  to  "thin" 
plates,  the  precision  remains  within  ±0.15X  which  Is  considered  excellent. 
Similar  results  were  obtained  with  an  epoxy  specimen  [5). 

Finally,  we  subjected  our  technique  to  a  very  difficult  test  of  measuring 
wavespeed  and  attenuation  In  a  highly  dispersive  as  well  as  attenuative  medium 
namely,  a  random  particulate  composite.  This  material  has  been  studied  In 
detail  by  Klnra  et  al  (7,8),  therefore,  we  had  a  database  with  which  we  could 
compare  our  present  measurements.  The  results  are  shown  in  Fig.  5.  Here  <Cj> 
Is  the  wavespeed  In  the  composite,  Cj  Is  the  wavespeed  In  the  neat  matrix 
material,  o  Is  a  suitably  normalized  frequency  and  c  Is  the  volume  fraction  of 
inclusions.  The  Important  point  Is  that  the  entire  dispersion  curve  (<cp 
vs  q)  as  well  as  the  attenuation  curve  was  obtained  In  a  single  experiment. 
Previously,  ( ( 7,8 J  for  example)  each  data  point  required  a  separate  toneburst 
experiment.  Therefore,  the  new  technique  Is  not  only  more  time-efficient  by 
several  orders  of  magnitude.  It  Is  also  more  accurate  by  about  one  order  of 
magnitude. 
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TABLE  1.  Test  Results  on  Aluminum  Sample 


Material: 
Wave  Type: 
Mode: 

Frequency: 

Oensity: 

Aluminum 

Longitudinal 

Transmission 

10  MHz 

2.8177  t  0.0004  g/ml 

h 

h/x 

C 

o/c 

Technique 

nm 

um/ysec 

% 

2.807 

4.4 

6.3572 

Tone burst 

2.807 

4.4 

6.3239 

0.013 

Second/First 

2.807 

4.4 

6.3275 

0.010 

All/First 

1.686 

2.7 

6.3461 

0.040 

All/First 

1.001 

1.6 

6.3538 

0.030 

All/First 

0.613 

0.96 

6.3594 

0.130 

All/First 

0.258 

0.4 

6.3231 

0.140 

All/FIrst 
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Fig.  5  Velocity  and  attenuation  versus  frequency  for  a  highly 
dispersive  medium. 


2.4  Throuqh-the-Thic<ness  Measurements.  Results 

We  now  present  some  results  concerning  the  assessment  of  damage  using  the 
techniques  developed  during  this  investigation.  The  attention  thus  far  has 
been  confined  to  cross-ply  laminates.  both  the  "thick'1  and  the  "thin" 
laminates  have  been  tested.  The  principal  conclusions  are: 

1.  As  expected,  there  is  no  measurable  change  in  the  through-the- 
thlckness  stiffness  due  to  transverse  cracking 

2.  The  attenuation,  on  the  other  hand.  Is  a  very  sensitive  and 

reliable  measure  of  transverse  cracks. 

Following  lay-ups  were  tested:  [  [0g,904,02ls,  and 

[0,9041s.  The  laminates  were  subjected  to  monotonic  tensile  loading  and 

transverse  cracks  developed  In  the  90°  plies.  The  "damage"  was  quantified  In 
the  form  of  edge  replication;  see  Fig.  6.  After  each  of  the  six  load  steps, 
the  coupon  was  removed  from  the  tensile  testing  machine  and  subjected  to  an 
ultrasonic  examination.  In  Fig.  6  we  have  plotted  the  change  in  attenuation 
from  its  undamaged  state,  Ak^x,  against  the  load  level;  the  damage  state  at 
each  load  level  is  included  schematically.  We  have  shown  only  that  one  Inch 
segment  of  the  length  of  the  coupon  which  is  Insonlfled  by  the  ultrasonic 

beam.  It  Is  well-known  that  the  scattering  of  a  wave  by  a  defect  Is  a 

sensitive  function  of  the  wavenumber,  k^a,  when  k^a  is  of  the  order  of  one, 
where  "a"  for  a  crack  is  the  half-crack  length.  To  study  the  influence  of  k^a 
we  conducted  tests  at  three  different  frequencies,  namely,  2.25,  5.0  and  7.5 
MHz  where  kja  *  1.25,  2.7  and  4.0,  respectively.  In  Fig.  6  we  observe  that 
Ak£X  Increases  dramatically  with  damage,  it  Is  therefore  a  sensitive  measure 
of  damage.  In  particular,  we  note  : nat  In  going  from  damage  state  3  to  4, 
only  one  additional  crack  appears.  the  Increase  in  attenuation  Is  quite 

significant  and  substantially  greater  than  the  error  of  measurement  which  is 
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Fig.  6  Attenuation  versus  damage  (transverse  cracking)  at 
one  location  and  three  frequencies. 
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estimated  to  be  about  t5<  (kgX  Is  of  the  order  of  0.2). 

The  damage  In  composite  materials  Is  not  distributed  uniformly;  In  fact, 
by  now  it  is  fairly  well  accepted  that  It  Is  a  stochastic  phenomenon.  In  view 
of  this  we  monitored  damage  at  three  different  locations  along  the  length  of 

the  beam  with  the  following  extremely  important  objective  In  mind.  Damage  is 

a  function  of  the  history  of  loading.  Yet  our  technique  should  measure  only 
the  current  state  of  damage  and  should  be  independent  of  the  past  history. 
Therefore,  In  Fig.  7,  we  have  plotted  attenuation  versus  damage  at  three 
different  locations  and  at  three  different  frequencies;  here  damage  Is  defined 
as  the  cumulative  crack  length  l.e.  number  of  transverse  cracks  times  the 
crack  length  which,  for  the  present  case  of  jOgSO^ls  Is  the  thickness  of  the 
904  group  of  plies  l.e.  0.5  mm.  Indeed,  the  three  curves  coincide  within  the 
errors  of  measurement.  Thus  the  curves  in  Fig.  7  may  be  viewed  as  Master 

Curves  for  an  Inversion  of  the  attenuation  data  to  yield  a  measure  of  damage. 

The  corresponding  data  for  the  wavespeed  at  one  position  Is  included  in 

Fig.  8.  There  Is  virtually  no  change  In  the  wavespeed  with  damage,  (this  Is 

consistent  with  a  well-known  fact  In  fracture  mechanics  that  a  crack  causes  a 

minimal  reduction  in  the  stiffness  of  a  structure  In  the  direction  parallel  to 

the  crack  faces).  Finally,  for  brevity  we  Include  here  only  the  final 

conclusions  of  some  related  Investigations.  The  details  may  be  found  in  Ref. 
(9). 

1.  It  has  been  reported  at  conferences  that  by  using  the  Stress  Wave 
Factor  approach  one  can  predict  the  site  of  the  final  fracture. 
8y  using  the  attenuation  as  a  damage  metric  we  found  the  converse 
to  be  true  l.e.  one  cannot  predict  the  site  of  the  final 

fracture.  We  hasten  to  add  that  our  observation  is  consistent 

with  the  theoretical  work  of  Riefsnider  et  al  at  VPI  &  SU  1101. 
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Fig.  8  Wavespeed  as  a  function  of  applied  stress  for  a 
[0/-.90- ,0o]„  laminate. 


2.  There  has  been  some  debate  In  the  literature  whether  or  not  the 
process  of  unloading,  subsequent  to  the  formation  of  transverse 
cracks,  can  cause  partial  (or  even  total)  crack  closure.  Our 
measurements  show  clearly  that  some  crack  closure  does  take  place 
191. 


2.5  Lamb  Wave  Method.  Theory 

An  extensive  literature  review  concerning  the  propagation  of  Lamb  waves 
In  an  anisotropic  plate  Immersed  In  a  fluid  was  carried  out.  Only  approximate 
solutions  were  found;  these  have  been  reviewed  in  [11].  For  the  case  in  hand 
the  density  of  water  is  comparable  to  the  density  of  the  composite  material. 
Therefore,  the  approximate  theories  -  which  are  based  on  the  assumption  that 
the  density  of  the  fluid  is  small  compared  to  the  density  of  the  plate  -  yield 
unacceptably  erroneous  results.  This  defined  the  first  objective  of  our 
work :  To  obtain  an  exact  analytical  solution. 

Consider  a  symmetric  composite  laminate  (not  necessarily  balanced).  Let 
xj  and  X£  He  In  the  plane  of  the  plate;  let  X3  be  perpendicular  to  It.  The 
plate  occupies  the  space  X3  *  tb.  The  generalized  Hooke's  law  for  such  a 
material  Is 
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This  relation  can  be  written  In  the  indlclal  form  as 
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°ij  =  c1jkl£kl 


for  1,j,k,l=1.2,3 


where  the  strain  displacement  relation  can  be  written  as 
e1j  3  (U1,j  +  UJ,1)/2  for  i.j-1.2.3 

The  equation  of  motion  In  an  elastic  medium  is 

°uj‘  oUi 

The  boundary  conditions  are: 

o31  *  0  and  o 33  *  -p  at  the  plate  faces  2  *  ±b. 

Here  p  Is  the  pressure  at  the  plate/water  interface.  The  wave  motion  in  water 
satisfies  the  equation  of  motion 


The  continuity  of  the  normal  displacements  requires  that  the  displacement  In 
solid  (U3)  be  equal  to  that  In  the  fluid  (WL): 

U3  =  WL  at  z  «  tb 

We  have  solved  this  boundary  value  problem  for  the  complex-valued 
wavenumber  kx  whose  real  part  (kj)  is  connected  to  the  phase  velocity  of  the 
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Lamb  wave  through  the  relation  c  =  u/k^  and  whose  imaginary  part  (k£)  is  the 
attenuation  of  the  Lamb  wave  l.e.  the  wave  decays  at  expf^x)  in  the  x 
direction.  We  point  out  that  this  is  not  attenuation  in  the  sense  of 
absorption,  for  both  media  are  elastic;  rather,  the  wave  attenuates  because  it 
leaks  energy  into  water  as  it  propagates.  It  Is  for  this  reason  that 
frequently  these  waves  are  called  "leaky  Lamb  waves."  We  now  present  the 
final  results  of  our  calculation.  As  with  Lamb  waves  for  a  plate  In  air,  we 
get  symmetrical  (longitudinal)  and  ant  1-synmetr leal  (flexural)  waves.  The 
wavenumber  kx  Is  a  root  of  the  transcendental  equations. 

I  M  B„  .  fc.1  „  A 

Tln(kzmb*  Gl"*Hp  0(Gm/o)SQR(k^-k^)  Hp 

Symmetric  Mode 

Tan(k  b)  Gp*Hm  1p,u2Cot(k  b)  .. 

t-7.-zw-k\  -  t - i2—  .  rp  _  r»]  »  0. 

Tan(k2pb)  6m*Hp  p(Gm/0)SQR(k^-kL)  Hp 

Asymmetric  Mode 

For  details  of  the  calculation  the  reader  Is  referred  to  1121.  The  roots 
of  these  equations  had  to  be  searched  numerically  and  that  was  found  to  be  a 
fairly  Involved  process.  Briefly,  the  modified  Newton's  (or  Secant)  method 
was  used  In  two  steps;  (1)  First,  assume  k  is  real  and  find  the  best  estimate 
of  k;  (2)  Second,  let  k  be  complex,  k  =  +  1kg.  Each  of  the  above  (complex) 

equations  yields  a  pair  of  real  equations.  A  simultaneous  numerical  solution 
Is  then  found. 
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2.6  Lamb  Wave  Method.  Experimental 

A  schematic  of  the  apparatus  Is  shown  In  Fig.  9.  In  theory  we  assumed 
time  harmonic  motion;  In  order  to  simulate  that  In  a  laboratory  we  use  a 
toneburst  l.e.  a  sinusoidal  wave  of  about  20  cycles  duration.  All 
measurements  are  made  away  from  the  two  ends  of  the  toneburst  (l.e. 
transients)  and  near  the  center  of  the  toneburst.  The  heart  of  the  system  Is 
a  pair  of  accurately  matched  transducers:  one  Is  energized  and  acts  as  a 
transmitter,  the  other  Is  passive  and  acts  as  a  receiver.  The  wave  travels  as 
a  longitudinal  or  pressure  wave  In  water.  The  specimen  is  inclined  to  the 
Incident  wave  at  the  correct  angle  for  the  Lamb  wave  mode  of  Interest  given  by 
the  Snel 1 ' s  law: 

Sin  1  S1n(»/2) 

Cw  " 

where  1  Is  the  angle  of  Incidence,  Cw  Is  the  speed  of  sound  In  water  and  CL  Is 
the  speed  of  the  particular  Lamb  wave  mode  of  interest.  The  leaky  Lamb  wave 
In  the  specimen  radiates  energy  into  water  syrnnetrically  on  both  sides;  this 
signal  Is  sensed  by  the  receiving  transducer.  The  attenuation  is  measured  as 
follows.  The  receiver  is  moved  gradually  in  a  direction  perpendicular  to  the 
"1 ine-of-sight"  of  the  transducers.  At  each  step  the  received  amplitude  Is 
recorded.  An  exponential  of  the  type  exp(-k2x)  Is  fitted  through  the  data  to 
obtain  k2-  The  velocity  is  measured  as  follows.  The  specimen  is  rotated 
gradually  and  In  very  fine  steps  and  the  correct  angle  of  Incidence  is  that 
which  produces  the  maximum  received  amplitude.  We  stumbled  upon  an 
Independent  (and  very  elegant)  cneck  on  the  correctness  of  the  angle  of 
Incidence.  Through  a  simple  calculation  jl2)  we  have  shown  that  if  one  moves 
the  receiver  perpendicular  to  the  "line  of -sight"  then,  for  the  correct  angle, 
the  arrival  time  of  a  point  of  constant  phase  at  the  receiver  remains 
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Fig.  9  Block  Diagram  of  the  experimental  setup. 


unchanged;  for  all  other  angles  the  arrival  time  changes,  the  magnitude  of  the 
change  being  proportional  to  the  error.  Experience  indicated  this  to  be  a 
sensitive  check  on  the  correctness  of  the  angle  of  incidence.  Having  found 


"1"  the  speed  Is  calculated  from: 
CL  =  Cy/Sin  (1) 


Calibration  Procedures  The  validity  of  our  theory  was  checked  against  some 
earlier  calculations  as  well  as  against  our  own  experimental  results.  Our 
theory  was  used  to  calculate  the  dispersion  (wavespeed  versus  frequency)  as 
well  as  the  attenuation  curves  for  a  steel  plate  In  water;  the  motivation  for 
choosing  steel  being  that  Merkulov  (131  had  earlier  reported  results  for  this 
material.  First  ten  modes  of  wave  propagation  were  calculated.  The  agreement 
with  the  theory  of  Merkulov  |13]  was  found  to  be  excellent;  see  (121  for 
details. 

Hext,  dispersion  and  attenuation  curves  were  calculated  for  an  aluminum 
plate  In  water;  wavespeed  and  attenuation  were  also  measured;  they  are 
compared  In  Fig.  10.  The  solid  lines  are  symmetrical  (longitudinal)  modes, 
the  dashed  lines  are  the  ani symmetrical  (flexural)  modes;  the  circles  are  the 
(discrete)  data  points;  the  agreement  between  theory  and  experiment  is 
considered  remarkably  good. 

Finally,  the  dispersion  and  attenuation  curves  for  an  undamaged  composite 
material,  (0,  903)5,  are  shown  in  Fig.  11.  It  was  found  that  It  was  easier  to 
collect  data  with  the  fundamental  symmetric  mode  ($0).  The  comparison  between 
the  theory  and  the  experiment  is  considered  quite  satisfactory.  In  the  next 
section  we  will  apply  the  Lamb  wave  method  to  the  assessment  of  damage. 
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2.7  Lamb  Wave  Method.  Results 

Three  types  of  laminates  were  studied:  [02*9020 1  s ,  [0,903)5  and 
[0.904]s.  The  form  of  damage  produced  was  transverse  cracking.  As  before, 
the  coupons  were  loaded  In  an  Instron  testing  machine  to  various  load 
levels.  At  each  load  step  the  test  was  Interrupted,  an  edge  replication  was 
taken  to  obtain  a  record  of  the  transverse  cracks  and  the  coupon  was  subjected 
to  an  ultrasonic  examination.  The  loading  was  then  resumed.  Several  modes  of 
Lamb  wave  propagation  were  examined.  Experience  Indicated  that  the  most 
suitable  mode  was  the  fundamental  symmetric  mode  SQ  (14).  All  data  presented 
In  the  sequal  were  gathered  by  using  this  mode.  It  also  turns  out  that  in  the 
long-wavelength  limit  the  dispersion  equation  for  this  mode  simplifies  to 

C2  =  Eu/p(1-v12v21) 

and  since  for  the  particular  lay-ups  under  consideration 

v12v21<<1* 

2 

approximately  C  =  E^/o 

Results  concerning  [02»902»0) s  are  reported  first,  see  Fig.  12.  Note  a 
rather  high  crack  density.  The  crack  occupies  the  90°-ply  group  l.e.  2a  = 
0.25  mn.  Because  the  crack  length  is  so  small  It  Is  not  surprising  that  the 
reduction  In  stiffness  Is  very  small,  about  5X  up  to  failure.  In  Fig.  12,  E0 
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Figure  12  Reduction  in  stiffness  in  a  [(^OgO] 
laminate,  f=0.5  MHz,  da0.55  mm. 
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Is  the  initial  stiffness  and  E  Is  the  “current"  stiffness  of  the  damaged 
specimen.  Next  we  report  results  obtained  with  [0,9031  s*  see  Fig.  13,;  the 
major  difference  here  Is  that  the  90°-group  has  six  plies  l.e.  2a=0.75  mm.  As 
a  result  a  relatively  greater  reduction  In  stiffness  Is  observed;  about  15%  to 
failure.  The  attenuation  measurements  are  also  reported;  attenuation  Is 
rather  small.  Next,  we  consider  a  { 0,90^} s  laminate;  see  Fig.  14.  Here  the 
edge  replications  are  also  Included  as  Is  the  location  of  the  transmitting 
transducer  (T)  and  the  receiving  transducer  (R).  Note  the  high  density  of 
transverse  cracks  at  the  (final)  load-step  9.  As  expected  on  the  basis  of  the 
larger  crack  length,  the  reduction  in  stiffness  is  very  large,  about  30%  to 
failure. 

Of  course  a  crack  may  be  viewed  as  "small"  or  "large"  depending  upon  its 
relation  to  the  wavelength.  To  this  end  we  have  defined  a  dimensionless 
wavenumber  k^a  =  2*fa/CL  where  f  Is  the  frequency  and  CL  Is  the  Lamb  wave 
speed.  The  data  collected  with  the  three  different  laminates  Is  "unified" 
Into  a  single  plot  in  Fig.  15.  The  cumulative  crack  length  Is  the  number  of 
cracks  times  the  crack  length,  2a.  It  Is  well-known  that  at  the  very  low 
wavenumber  of  k^a  =  0.06,  we  are  In  the  so-called  Rayleigh  regime  of 
scattering  where  the  cracks  act  as  weak  scatters.  As  a  result  attenuation  Is 
very  low.  At  a  slightly  higher  K^a  =  0.28,  the  attenuation  Is  still  rather 
low.  But  when  kja  *  0.45,  we  get  a  substantial  attenuation.  It  Is  also  well- 
known  that  at  certain  critical  frequencies  cracks  begin  to  resonate  (very  much 
like  the  vibrations  of  a  string)  This  typically  happens  when  k^a  =  0(1). 
Figure  15  suggests  that  we  may  be  approaching  resonance.  This  conjecture  Is 
the  subject  of  the  current  Investigation. 

In  conclusion,  we  have  established  that  both  the  wavespeed  and  the 
attenuation  of  Lamb  waves  serve  as  measures  of  damage  in  composite  material. 
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gure  13  (a)  Reduction  in  stiffness,  (b)  increase  in 

attenuation  in  a  [0  903]s  laminate,  f«0.5  MHz, 
ds0.55  mm. 


CUMULATIVE  CRACK  LENGTH  (mm) 

Figure  is  Increase  In  attenuation  (((2)  with  cumulati¬ 
ve  crack  length  for  three  laminates  tested. 
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